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Abstract 

We consider the two-loop self-mass sunrise amplitude with two equal masses M and the external 
invariant equal to the square of the third mass m in the usual d-continuous dimensional regularization. 
We write a second order differential equation for the amplitude in x — m/M and show as solve it in close 
analytic form. As a result, all the coefficients of the Laurent expansion in (d — 4) of the amplitude are 
expressed in terms of harmonic poly logarithms of argument x and increasing weight. As a by product, 
we give the explicit analytic expressions of the value of the amplitude at x — 1, corresponding to the 
on-mass-shell sunrise amplitude in the equal mass case, up to the (d — 4) term included. 



1 Introduction 

It is known that the sunrise graph of Figure |l| in the arbitrary mass case has four master integrals; they 




Figure 1: The sunrise graph for arbitrary masses. 



can be taken to be the fully scalar amplitude 
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where d is the continuous dimension, irii,i = 1,3 the masses, p 2 the square of the external momentum (all 
momenta being Euclidean), and the three amplitudes Fi defined as 

Fi{d,m\,m\,m\,-p 2 ) = --^- s F (d, m\, ml, m\,p 2 ) , (1.2) 

for i = 1,3. The higher mass derivatives of the amplitudes can all be expressed, as a consequence of the 
integration by parts identities, in terms of the four master integrals. The four master integrals satisfy 
further a system of linear non homogeneous differential equations in j> 2 , derived in Q. 

In the particular kinematical configuration — p 2 = m\ = m 2 and m\ = m 2 = M 2 , Figure 2, one has 




Figure 2: The sunrise graph with two equal masses and the external invariant equal to the third squared 
mass. 



F (d,m 2 ,M 2 ,M 2 ,p 2 = - m 2 ) = J ^ 2 _ 2 _ _ _____ _ _L _ __ _ _ ■ (1 . 3) 
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introducing the dimensionless variables x = m/M one can define the dimensionless function $(d, x) as 

F (d,m 2 ,M 2 ,M 2 ,~m 2 ) = M 2d - 6 C 2 (d)<S>(d, x) , (1.4) 

i.e. 

$(d,x) = C- 2 (d)M 6 - 2d F (d,M 2 x 2 ,M 2 ,M 2 ,-M 2 x 2 ) , (1.5) 

4— d 

where C(d) = (47r)~2~r(3 — d/2) is an overall normalization factor, with the limiting value C(4) = 1 at 
d = A. 

It will be shown in this paper that the integration by part identities for the master integrals imply 
a second order differential equation in x for x); when functions and equations are systematically 
expanded around d = 4 as a Laurent series in {d — 4), one obtains a system of chained differential equations 
which can be easily solved recursively in closed analytical form, virtually up to any order in (d — 4); one 
finds in this way that the coefficients of the (d — 4) expansion of $(<i, x) are a combination of harmonic 
polylogarithms || of increasing weight. 

The differential equations themselves with the knowledge of the behaviour at m = can be used to fix 
completely the integration constants; as a byproduct we give the value of the integral at m — M (on mass 
shell equal mass case) up to the order (d — 4) 5 included. 

2 The differential equations. 

For deriving the differential equation for <5>(<i, x) we introduce the auxiliary functions 

^ Q (d,r],p 2 ) = m 6 - 2d F (d,m 2 ,m 2 r],m 2 T],p 2 ) , (2.1) 
*i(d,?7,p 2 ) = m s ' 2d [F2(d,m 2 ,m 2 ri,m 2 ri,p 2 ) + F 3 (d,m 2 ,m 2 r 1 ,m 2 ri,p 2 )] . (2.2) 

One finds at once 

1^0(^7?,/) = 77, P 2 ) , (2.3) 

from which in turn 

^-* 1 (d, ?? ,p 2 ) = -m( 10 - 2d ) 
or] 

d \ 

+ g^|J [F 2 (d,m 2 ,m 2 ,mlp 2 )+F 3 (d,m 2 ,mlmlp 2 )] ml=m2 _ m2ri .(2.4) 
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As already remarked, the mass derivatives of the master integrals Fi,i — 1,3 can be expressed, by using the 
integration by parts identities, in terms of the master integrals themselves (see for instance the Appendix 
of |ffl)), so that the above quantity is equal to a combination of F , corresponding to $ (d, rj,p 2 ), of the sum 
F 2 + F 3 (as required by the m 2 = symmetry), corresponding to \&i (d, 77, p 2 ), and of the fourth master 
integral F 1; independent of 'J'o and P ms an inhomogeneous part given by a few products of tadpoles 
T(d,m 2 ),T(d,M 2 ), with 



m -') - / (2 ')L k ,l m , = CM (<i 4) ■ ( "> 

The coefficients of the various terms are ratios of polynomials in d, p 2 and the masses; in the limit p 2 —> —m 2 
the denominators of the coefficients develop a singular l/(p 2 + m 2 ) 2 behaviour, which is however fully 
compensated by corresponding zeroes in the numerators, so that the limit p 2 — —m 2 is finite. In that 
limit, further, the coefficient of the master amplitude F\ vanishes, so that -j^'fy \{d,ri,p 2 — — m 2 ), which is 
essentially the second 77-derivative of ^>o(d,r],p 2 = —m 2 ), is expressed in terms of ^i(d, rj,p 2 = — m 2 ), its 
first ?7-derivative (up to a sign), and of ^o{d, ?7,p 2 = — rn 2 ) itself, plus an inhomogeneous term given by the 
tadpoles. For taking the p 2 = —m 2 limit, let us introduce the function ^(d, rj) through 

F (d,m 2 ,M 2 ,M 2 ,-m 2 ) = m 2d ~ 6 C 2 (d) *(d,»/) , (2.6) 
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i.e. 



^{d,r}) = C~ 2 (d) m & - 2d F (d,m 2 ,m\m\ -m 2 ) . 
In terms of ^(d,rj) the second order differential equation for *Sf (d,r],p 2 — —m 2 ) then reads 



(2.7) 



d 2 

V(l-V)^(d,r])- 



|-|» 7 +(l-2» ? )(d-4) 



l + \{d-A) + - A {d-A) 2 



ry) 



v (d-A)/2 7] l+(d-4) 



4(d-4) 2 8(d-4) 2 



(2.8) 



The equation, exact in the continuous dimension d, has been written in terms of powers of (d — 4) for 
convenience of later use. 

As r\ = l/x 2 , by comparing the definitions Eq.s (fOL |2.6|) one has 



* [d,Ti 



1 



x b ~ 2<1 $(d,x) , 



(2.9) 



and the above equation for ^(d, rf) can be converted in the following equation for <I?(d, x) 



~(l-z 2 ) + [\ -\xA (d-4) + i(d-4) 2 



$(d,a:) 



x 



2+(d-4) 



4(rf - 4) 2 8(d-4) 5 



(2.10) 



3 The expansion in (d — 4) around d = 4 . 

It is known [fij that in the d — ► 4 limit $(d, cc) develops a double pole in (d— 4), so that its Laurent expansion 
in (d — 4) reads 



$(d,x) = ]T (d-4)"$^(a 



(3.1) 



By inserting the above expansion in Eq. ( 2.10| ) one obtains a system of chained inhomogeneous equations 
in the functions & n '(x). For n — —2, —1 they read 



d 



1 



1 



d__2_ 

dx 2 \x 1 — x 1 + x J dx x 2 
d 2 (2 1 1 \ d 2 



dx 2 \x 1 — x 1 + x J dx x 2 
3 1 1 



x 2 2(1-2;) 2(1+ x) 



$(- 2 )(x) 



1 



1 



2x 2 4(1 -x) 4(1 + a:) 
1 



1 — a; 1 + s / 9i 
1/1 1 



' ^* ( - 2) (x) 



2 VI -x 1 + x 



Inx 



(3.2) 



while for n > one has 
^ /2 1 



drz; 2 



x 1 — x 1 + x J dx x 



$(«)(V) 



1 



d 



+ {x 2 + 2(1 - x) + 2(1 + x)) ^ 1){X) + {x 2 



1 — x 1 + x J dx 
1 1 1 



x 2 2(1 -x) 2(1 + x) 
1(1 1 \ ln n+2 a: 



$(«-!) (a;) 

$("~ 2 )(x) 



2 \l-z ' 1 + ^7 (rc + 2)! 
A few comments are in order. Quite in general, all the equations have the form 

dx 2 + \x + l-x 1 + x dx x 2 [) H W ' 



(3.3) 



(3.4) 
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the homogeneous part of the equations is always 



dx 2 



1-x 1 



d_ 

dx 



ct>{x) = , 



(3.5) 



while the inhomogeneous terms (x) vary from equation to equa tion and involve, besides the coefficients 
of the (d — 4)-expansion of the inhomogeneous term of Eq.( [2.ic| ) which are explicitly known since the 
beginning, also the coefficients ^ n ~ 1 \x) and $(™~ 2 )(:e) of the (d — 4)-expansion of the unknown function 
^(d, x). However, in a systematic bottom- up approach to the solution starting from the smallest value of 
n, i.e. n = —2, the inhomogeneous part can also be considered as known. 

Finally, all the coefficients appearing in the equation, both in the homogeneous and the inhomogeneous 
part, are simple combinations of the rational factors 1/x, 1/(1 — x) and 1/(1 + x); powers of \nx are also 
present, from the expansion of the original inhomogeneous term. 



4 The behaviour of the equation for x — > 



+ 



As x is proportional to a mass and therefore is an essentially real and positive variable, we study Eq.( [2.1C ) 
in the x — > + limit; to that aim, we try the most general expansion 



<!>{d,x) 



\n=0 



(d)x r ' 



(4.1) 



By substituting the above expansion in Eq.( 2.10 ), we find for the leading exponents oti the four allowed 
values ql\ = (2 — d), cti = (d — 3), corresponding to the two independent solutions of the associated 
homogeneous eq uatio n and az = (d — 2), a± = 0, as required by the inhomogeneous terms appearing in 
the r.h.s. of Eq. (|2.10|) . 

By direct inspection of Eq. (|l.3|) and Eq.(L4), one sees that $(d, x) is finite for d > 2 in the x — > + limit, 
so that the first two behaviours, with exponents (2— d) and (d— 3), are ruled out (i.e. the coefficients a« (d) 
are all vanishing for i = 1,2 and any value of n ). The coefficients of the other terms are then completely 
fixed by the inhomogeneous part of the equation, and an elementary calculation gives the following explicit 
behaviour in the x — > + limit 



1 



2(d-2)(d-4) 2 



1 



2x 2 



d-3 d(d-5) 



+ x 2 x d - A 



0(a 



(4.2) 



When expanding in (d — 4) according to Eq.(p4|), the above equation reads 

1 1 



$(-2)( 
$(-!)( 



Oix* 

,,.2 



16 



-lna; + 0(x 3 



& \x) = -i + ^x 2 + ^-(Inx - ln 2 x) + 0(x 3 ) , 



(4.3) 



and so on up to any desired order in (d — 4) . 

It is to be observed that by using the differential equation Eq.( 2.10| ) one can transform a qualitative 
information, the finiteness of $(<i, x) at x = 0, into a quantitative information, such as the actual value of 
the x — > + behaviour of $(d, x) given by Eq.(4.2), without carrying out any explicit loop integration in 
Eq.O). 
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5 The solution. 



The algorithm for solving Eq.(3.4) go es b ack to Eulcr. Let </>i(x), 4>i{x) be two independent solutions of 
the associated homogeneous equation (3J3), such that their Wronskian W(x), defined as 



w(x) = M^ix) - 4>[(x)Mx) , 



(5.1) 



does not vanish; the solution of the second order Eq^^) can then be written as 

where , $2 are two integration constants, to be fixed by the boundary conditions Eq.(4.2) expanded 
in (d-A), i.e. by Eq.sQ. 

For the actual use of Euler's algorithm, the solutions of the homogeneous equation (3.5) are needed, ft 
turns out that the homogeneous equation associated to Eq. (|Tq) at d — 4 



r)(l - r)) 



dr\ 2 



d_ 

di] 



(5.3) 



is a hypergeometric equation, whose solution regular at rj — is the very simple hypergeometric function 



2Ft 



1 „ 1 



(5.4) 



With that input and the transformation (2.9) it is easy to obtain the following solutions of the homogeneous 
equation (EO): 



having the Wronskian 



and the 



behaviour 



4>i{x) 
(j> 2 {x) 



(1 - X 2 ) 2 
x 2 

(l-x 2 ) 2 



hi 



W(x) = 16- 



1 -x 



(5.5) 



(5.6) 



cj) 2 {x) 



-6x 



(5.7) 



With the above explicit homogeneous solutions one can solve recursively the system of the non homogeneous 
equations (3.3) for &( n \x), starting from n — —2, by carrying out repeatedly the integration in x' of (5.2) 
in close analytic form; one obtains in that way a combination of Harmonic PolyLogarithms (HPLs) of 
increasing weight [Q times simple coefficients involving the rational factors 1/x, 1/(1 — a;) and l/(l + x). 

To illustrate this statement, which is the central point of this paper, let us recall for convenience of the 
reader the definition of the HPL's. The HPL's form a family of functions depending on an argument, say 
x, and on a set of indices, say ai,i = 1, .., w or a in more compact notation, where each of the can take 
one of the three values 1,0,-1 and whose number w is called the weight of the HPL. At weight w = 1 there 
are 3 HPL's, defined as 

H(l;x) = -ln(l-x), 
H(0;x) = lnx , 
H(-1;.t) = ln(l+ar) , (5.8) 



5 



whose derivatives can be written as 

— H(a;ar) = f(o,x) , (5.9) 
ax 

where the index a can take one of the three values (1, 0, —1) and the rational factors f(a; x) are given by 

1 



f(l;z) 
f(0;x) 



x 

f(-l;x) = jl-. (5.10) 

At weight w > 1, if all the u> indices are equal to let us indicate them by W and define correspondingly 

H(6 w ;x) = ^-\n w x; (5.11) 

in all the other cases (i.e. when the indices are not all equal to zero), let us indicate any set of w indices 
by (a, b), where a can take one of the three values (1, 0, —1) and b stands for the set of other w — 1 indices, 
and define correspondingly 

R(a,b;x)= [ dx'f(a; x')U(b; x) . (5.12) 
Jo 

Note that in full generality one has 

4-H(o, 6; x) = f(o, x)R(b; x) , (5.13) 
ax 

which can also be written as the equivalent indefinite integration formula 

J dx'{(a;x')R(b;x') = A + R(a,b;x) , (5.14) 

where A is an integration constant. 

Further, the product of two HPL's of a same argument x and weights p, q can be expressed as a 
combination of HPL's of argument x and weight r = p + q, according to the product identity 

K(p;xMq;x) = £ K(r;x) , (5.15) 

r—p^±lq 

where p,q stand for the p and q components of the indices of the two HPL's, while pttl q represents all 
possible mergers of p and q into the vector r with r components, in which the relative orders of the elements 
of p and q are preserved. The simplest cases of the above identities are 

H(a; x)H(b; x) — H(a, b; x) + H(6, a; x) , 
K(a; x)H(b, c; x) = H(a, 6, c; x) + H(6, a, c; ie) + H(&, c, a; ; (5.16) 

more complicated cases are immediately established recursively (all the above formulae can indeed easily 
be checked by differentiating, repeatedly when needed, with respect to x). 

After this digression, let us come back to Eq.(5.2), which gives the solution as an indefinite integral. 
The original n = — 2 integrand consists of powers of a;', 1/(1 — x') and 1/(1 + x'), (without any loss of 
generality (1 — x') and (1 + x') can be taken to occur only as negative powers), which are present in all 
the factors (see the first of Eq.s([T^) for the definition of i?( -2 ) (x 1 )), and of the two HPL's of weight 1, 
H(l; x') = — ln(l — x') and H(— 1; x') — ln(l + x'), which occur in 4> 2 {x'), Eq.( |5.5|) . The terms containing 
powers of the rational factors different from the mere inverse, i.e. different from the first powers of 1/x', 
1/(1 — x') and 1/(1 + x'), can be integrated by parts; in so doing they give rise at most to derivatives of 
the occurring HPL's, i.e. again products of the three rational factors x', 1/(1 — x') and 1/(1 + x'), and, 



G 



according to Eq.s fl5.13j 5.9), of HPL's of lower weight if the HPL's in the integrand have weight greater 
than 1 (as it is in the general case n > —2). The result can be partial fractioned, and the integration by 
parts repeated until all powers of the rational fractions different from the mere first power are processed. 
The terms with the mere inverse of the three rational fractions and HPL's, finally, can be integrated at 
once according to Eq.( 5.14 ), giving ri se to HPL's of higher weight. The integration constants are then fixed 
by imposing the x — > + behaviour (4.2). It is clear that the argument applies to any of the subsequent 
integrations as well, recursively in n. 



(—2) ( — 2) 

For n — —2 the explicit calculation gives §\ = 0, § 2 =0 and 



( - 2) (*) = 



(5.17) 



Note that, in this case, all the HPL's cancel out from the final result. 

One can then proceed to n — — 1; the term lnx appearing in the inhomogeneous part of the second 
of Eq.s(3 : 2) can be written as H(0;x), see (|5.8[ ), the products of the HPL's already present in i?( _1 )(x') 
an d the HPL's of the <j>i{x') can be rewritten as combinations of single HPL's of suitable indices according 
to 5.15, the integration performed as explained above - and so on for the higher values of n, up to any 



required order. An explicit calculation gives for the integration constants &\ 



(-i) 



(i) 



107/1024, $f ) = 185/4096, ®\ 



(3) 



while, for any n, <&2 = 0- The first <I>( n )(:z;) are 



-1285/16384, 



(4) 



-5/64, $< 0) = 29/256, 
18991/65536, <^ 5) = 164173/262144, 



(5.18) 



(l-x 2 ) 2 



±L + ^H(0;x)- ±H(0;x)H(Q;x) 



8x 2 



H(0;x)H(-l]x)-H(0;x)H(l;x)-H(0,-l;x) + H(0,l;x)) , (5.19) 



(l-x 2 ) 2 

J; v 

32 x 2 

( 5 [ H(0; x)H{-l; x) - H(0; x)H(l; x) - H(0, -1; x) + H(0, 1; x) 

H(0; x)H(0; x)H(l;x) - H(0; x)H(0; x)H(-l; x) - H(0; x)H(-l; x)H(-l; x) 
+ 4iT(0; x)H(-l; x)H(l;x) - H(0; x)H(l; x)H(l; x) 
4 [ H(0, -1; x)H(-l; x) + 2H(0, -1; x)H(0; x) - 2H(0, -1; x)H(l; x) 
4 [ H(0, 1; x)H(l; x) - 2H(0, 1; x)H(-l; x) - 2H(0, 1; x)H(0; x) 

H(0, -1, 1; a;) + H(Q, 1, -1; x) - |i?(0, 0, -1; a;) + ^iJ(0, 0, 1; x) 
-^(0,-1,-1;^) -±H(0,l,l;x) 



x 

(1 -xf 



2 r 



H(0, -1; x) - H{0; x)H(-l; x) 



H(0,l;x)-H(0;x)H(l;x) 



x 



55 11 

64 + ¥ 



H(0;x) 



15 13 
64 + 32 



if(0;x) 



16 



- -if (0; x)H(Q; x) + ^if (0; x)H(0; x)H(Q; x) 
H(0;x)H(0;x) . 



(5.20) 
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Eq.s( 5.17 , 5.18 ) are of course in agr eeme nt with the general mass case (see for instance 0, Q); Eq.( 5.1S ) 
reproduces Eq.(25) of j|, while Eq.( 5.2C ) is new in the literature. As matter of fact we evaluated explicitly, 
by means of an integration routine written in FORM ||, all the <&( n \x) up to n — 5 included, which is found 
to involve HPL's of weight 7. The explicit expressions obtained are however too long to be listed here. It 
is found that, in general, $ ( -™- ) (x) involves HPL's of weight (n + 2). 

The values of the HPL's of argument equal to 1 are known ||; we can then take our explicit analytic 
results for the & n ^(x) up to n = 5 and eva luate them at x — 1 by using ||, obtaining the analytical on- 
shell value of the sunrise amplitude Eq. 1.2 in the equal mass case limit up to the 5th order in the (d — 4) 
expansion. At x = 1 the terms of highest weight disappear from $w(a;), due to the general structure of 
the solution, Eq. ( |5.5| , j5~^ ), and the coefficient of (d — 4)™ involves constants up to weight n + 1 only. The 
result reads 



®(d,x = 1) 



17 1 



59 



(d-4) 2 32 (d- 4) 128 



+ (d-4) 



65 
512 



24 



id- 
id- 
id- 



4) 2 
4)3 

4) 4 



/ 1117 13 
V2048 ~ 96' 
13783 115 

384 
865 



-(d-4) 5 



1 



7r 2 ln2 



1_ 

16 



C(3) 



8192 
114181 



13 
32 
115 



_ — ^ 2 ln2 + — C(3) + -7r 2 ln 2 2 - 
'"• 1 64 



31 



-7T 



32768 1536 



■■ ■ - 2, „ 805 ,, n , 13 2 , ., 

7r 2 ln2 C(3) 7r 2 ln 2 2 

128 256 v ; 32 



2880 
2. 403 



11520 



1 i 4 9 3 
16 ln2+ 2 a4 

H ln 4 2 

64 



39 31 4l n 1 2l 3o 5 2 _. 3, 5n 9 465 A/r . 

a 4 7r 4 hi2 + -7r 2 ln 3 2 H tt 2 C(3) H ln 5 2 a 5 -{ C(5) 

8 960 8 96 sw 80 2 5 128 sv ; 



820495 5971 



131072 6144 



865 2l „ 
512^ 1112 



6055 115 2 2 

1024 C(3) + 128" ln 2 



345 

"32 
9 

160' 



■04 



65 
384 



^ 2 C(3) 



13 



403 



7r 4 ln2 



39 



ln & 2 



713 
9216 
6045 



H^ 2 

256 



C(5) 



117 



7r 4 ln 2 2 



1 

64 



7r 2 ln 4 2 



7T 2 ln 3 2 , . _ 

32 3840 320 512 

6 + i^C(3)ln 3 2 - !c(3h 2 hi2 + ^C 2 (3) 



79 



34560 



32 



128 



45 45 21 

-_C(5)ln2 - — a 5 ln2 + ln 6 2 

4 SW 4 5 320 



9a fi 



45 h 



+0((d-4) 6 ) , 



(5.21) 



where we have expressed the constant s6 appearing in jsj as 
s6 = 



ivr 2 C(3)ln2 - — 7r 2 ln 4 2 - — 7r 4 ln 2 2 + — tt 6 + 2C(5)ln2 
6 ^ W 72 720 540 sv ; 



-2a 5 ln2 + ~£(3) ln 3 2 + ^ln 6 2 + ^C 2 (3) - 4a 6 + 2b 6 



(5.22) 



and the definitions of all the previous constants in terms of HPL's or Nielsen's polylogarithms (NP1) is 
given in Table [l], t ogeth er with their numerical values. 

Our result, Eq. (|5.2l| ), agrees with previous results in the literature (apart from obvious normalization 
factors), such as [g), whic h goes up the (d — 4) 3 of present work, and which arrives at (d — 4) 4 . 

Numerically, Eq.(5.21) gives 



$(d, i=l) = -0.37500000000000000000- 



1 



0.53125000000000000000- 



'(d_4) a ' (d-4) 
-0.46093750000000000000+ 0.53818664171204166667(d- 4) 
-0.46186261021407291667(d- 4) 2 + 0.53810855601624843750(d- 4) 3 
-0.46190033063946289062(d- 4) 4 + 0.53809670843016515942(d- 4) 5 
+C((d-4) 6 ), (5.23) 



8 



in perfect agreement with the numerical results of which moreover gives two further terms in the 
expansion. 

As a last curiosity, let us observe that for d — > 4 Eq. ( 5.25 ) can be rewritten as 

$(d, x = 1) = , - 1 r ( +0.37500000000000000000 - 1 xo - 0.90625000000000000000—^— 

(d- 5) \ (d- 4) 2 (d-4) 

+0.99218750000000000000- 0.99912414171204166667(d- 4) 

+ 1.0000492519261145833(d- 4) 2 - 0.99997116623032135417(d- 4) 3 

+ 1.0000088866557113281(d- 4) 4 - 0.99999703906962805005(d- 4) 5 

+ 0((d-4) 6 )) , (5.24) 
showing that the coefficients of the powers of (d — 4) n approach quickly (— l) n for increasing n. 



constant 


HPL 


NP1 


numerical value 


C(3) 


ff(0,0,l;l) 




1.2020569031595942854 


0,4 


H(Q, 0,0,1; 1/2) 


5 3 ,i(l/2) 


0.51747906167389938633 


C(5) 


ff(0,0,0,0,l;l) 




1.0369277551433699263 


a 5 


#(0,0,0,0,1; 1/2) 


54,1(1/2) 


0.50840057924226870746 


a 6 


#(0,0, 0,0, 0,1; 1/2) 


ft,i(l/2) 


0.50409539780398855069 


h 


£T(0, 0,0, 0,1,1; 1/2) 


^4,2(1/2) 


0.0087230030575968884272 



Table 1: Constants up to weight 6 appearing in the calculation 



6 Conclusions. 

We have shown that the two-loop self-mass sunrise amplitude, in the usual d-continuous dimensional regu- 
larization and particular kinematical configuration corresponding to two equal masses M and the external 
invariant equal to the square of the third mass m, satisfies a second order differential equation in the 
dimcnsionless ratio x — m/M. The equation can be expanded in Laurent series at d = 4, and solved in 
closed analytical form. As a results, all the n-th coefficient of the Laurent expansion can be expressed in 
terms of harmonic poly logarithms of argument x and weight up to n + 2. As a by product, we give the 
explicit analytic expressions of the value of the amplitude at x = 1, corresponding to the on-mass-shell 
sunrise amplitude in the equal mass case, up to the (d — 4) 5 term included. 
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